Non-perturbative fixed point in a non-equilibrium phase transition 
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We apply the non-perturbative renormalization group method to a class of out-of-equilibrium 
phase transitions (usually called "parity conserving" or, more properly, "generalized voter" class) 
which is out of the reach of perturbative approaches. We show the existence of a genuinely non- 
perturbative fixed point, i.e. a critical point which does not seem to be Gaussian in any dimension. 
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Our understanding of equilibrium phase transitions is 
largely due to the success of perturbative renormalization 
group (RG) methods performed around some (upper or 
lower) critical dimension d c and to the existence of inte- 
grability and conformal symmetry properties in d = 2 . 
The situation is far less satisfactory out-of-equilibrium, 
where the relevant ingredients determining universality 
classes are sometimes not even known 0, ■ There are 
a number of technical reasons for this: (i) many systems 
possess neither a lower d c nor a low-dimensional exact 
solution; (ii) even at the critical point, models with a 
Langevin-like dynamics, that is, those that involve only 
one time derivative in their kinetic terms cannot be con- 
formal invariant; (iii) contrary to equilibrium, no RG cal- 
culation is in general available at and above three loop 
order and this prevents from using series re-summation 
techniques to compute accurately universal quantities in 
low dimensions. If, moreover, one keeps in mind that 
features not accessible to perturbative RG methods may 
play a crucial role, then the so-called non-perturbative 
renormalization group (NPRG) approach appears as a 
method of choice out-of-equilibrium. 

Such non-perturbative effects were evidenced recently 
in a study of the classic reaction-diffusion problem where 
particles A diffuse, branch (A — > 2A) and annihilate 
(2 A — > 0) with rates D, a, and A 0. Whereas, in an im- 
portant work [f| , Cardy and Tauber had shown that per- 
turbative RG calculations led to conclude that no finitc- 
<j transition to the empty absorbing state is possible for 
d > 2, an NPRG study at the non-universal level showed 
that such absorbing phase transitions exist in any finite 
dimension (with their critical properties in the directed 
percolation class, as expected from perturbative RG). 

In this Letter, we apply the NPRG method to a sim- 
ilar class of absorbing phase transitions, one for which 
even the universal properties are out of the reach of per- 
turbative approaches. We put forward the existence of 
a genuinely non-perturbative fixed point, i.e. a critical 
point which does not seem to be Gaussian in any di- 



mension. Our calculations unveil the structure of the 
RG flow, reveal clearly why perturbative methods are 
doomed to failure in this case, and provide estimates of 
critical exponents heretofore accessible only via numeri- 
cal simulations or series expansions. 

We actually treat two classes of systems known to be 
equivalent in d — 1, the physical dimension where the 
non-perturbative fixed point alluded to above is relevant. 
The first group 0,0 includes the reaction-diffusion sys- 
tem described above but where the branching reaction 
now creates pairs of particles (A — > 3A), so that, in- 
cidently, the parity of the total number of particles is 
conserved. Improperly named "parity-conserving" (PC) 
class (as argued in |8|, where this conservation law was 
shown to have no influence on similar reaction-diffusion 
systems) , it is best characterized by the second group of 
problems || , that of phase transitions into one out of two 
^-symmetric absorbing states ^(j- In d = 1, the par- 
ticles of the PC model can indeed be seen as interfaces 
between "+" and "— " domains 11]. In this "spin" lan- 
guage, domains evolve and compete under .^-symmetric 
rules with noise only acting at interfaces, the definition 
of the "generalized voter" (GV) class as given in [l^. 

The reaction-diffusion problem (A — > 3^4, 2A — ► 0) was 
also studied in |5|. Using the Doi-Peliti formalism [R 
one can obtain the following microscopic action 0, f 



S[tf>, 



0(d t -DV 2 )0-\{l~4> 2 )0 2 + a{l-4> 2 )W (i) 



in terms of the "physical" density field <p and the associ- 
ated response field <f>. Cardy and Tauber first performed 
an expansion around the upper critical dimension d c = 2 
where the transition occurs at zero branching rate a, so 
that the fixed point then is that of the pure annihilation 
problem 2 A — > 0. They showed that this annihilation 
fixed point remains relevant down to d = | (at one-loop 
order), where it becomes also attractive in the a direc- 
tion. In another expansion, this time performed directly 
in d = 1 , they were able to identify an appropriate combi- 
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nation of the coupling constants A, a which does admit a 
fixed point for d < |, although the flow diagram with re- 
spect to these original variables is rather peculiar. Even 
if their results do suggest that a phase transition should 
exist at a ^ — as observed in numerical simulations and 
mean-field-based expansions of many microscopic models 
in d = 1 0>@ — j the critical exponents remain poorly de- 
termined and, worse, the very possibility of computations 
beyond one loop-order appears to be problematic 

As for the GV class, the following Langevin equation 
was recently proposed [Tlj : 



d t <f)= (-a^ + ^ 3 )(l-<j) 2 ) + D\7 2 ( j)+y/2\(l -0 2 ) ry (2) 

with <fi E [—1,1] and rj is a delta-correlated Gaussian 
noise of unit variance. The (1 — <p 2 ) factors, appearing 
both in the deterministic force and in the noise ampli- 
tude, impose (j) = ± 1 to be symmetric absorbing states. 
In d = 1, only one type of transition is observed by vary- 
ing a for any value of /j,, and its critical properties are 
indeed those observed for the so-called parity-conserving 
models Q. Taking \x = for the sake of simplicity, 
the generating functional associated with the simplified 
Langevin equation is nothing but Q where <fi «-> cp and 
t — > —t. Thus, at this bare level, the two problems are 
strictly equivalent in d = 1. Wc shall see later that this 
equivalence is preserved under the RG flow. 

We cannot detail here the implementation of the 
NPRG but only mention the essential features of the 
method fl5t flq . The main idea is to build a one- 
parameter family of models, indexed by a momen- 
tum scale fc, interpolating smoothly between the short- 
distance physics at the (microscopic) scale k = A, where 
no fluctuation has been taken into account, and the long- 
distance physics at scale k — 0, where all fluctuations 
have been integrated out. In Wilson's original formula- 
tion, this leads to a flow of effective Hamiltonians — for 
the slow modes — defined at scale k. Following 0], we 
focus on the flow of "free energies" for the rapid modes, 
\q\ E [k, A], i.e. those already integrated out at this scale. 
This is achieved by adding a mass-like term of order k 2 
to the slow modes (\q\ < k), which "freezes" them. This 
mass term reads: 



AS fc [0,0]=/ R k tfM-q,-u)<t>(q,u) (3) 

J q,to 

where a convenient choice ^3 °f "cut-off" function is 

R k (q 2 ) = k 2 (l- q 2 /k 2 )e(l-q 2 /k 2 ) . 

The "partition functions" Z k [J,J] = J DcpD(pexp(-S - 
AS k + J J<f> + J J<p) become therefore /c-dependent. 
Through the Legendre transform of log Z k [J, J] , one ob- 
tains the state function T k — analogous to the Gibbs 
free energy at equilibrium — which depends on the fields 
tp = SlogZ k /SJ and $ = SlogZ k /SJ: 

r& [i>> $] + log z k [ J, J] = [ [ Ji>- [ R h W- (4) 



Note that the last term in Eq.Q, proportional to R k , 
ensures that T k has the proper limit at k — A: Ffc = A ~ S 
flij | . The following exact functional differential equation 
governs the RG flow of T k un der an infinitesimal change 
of the scale s = log(fc/A) [HEI: 



d s T k = -Tr 
2 



sR k (ri 2) 



Rk 



(5) 



where R k is the symmetric, off-diagonal, 2x2 matrix 

" (2) 

of element R k and 1^ [if), tp] the 2x2 matrix of second 
derivatives of T k w.r.t. ip and ip. Obviously, Eq. JSJ) can- 
not be solved exactly and one usually truncates it. A 
standard truncation is the derivative expansion |l5l | in 
which T k is expanded as a power series in V and d t . The 
local potential approximation (LPA), which is the sim- 
plest such truncation, consists in keeping only a potential 
term in T k while neglecting any field renormalization: 



r LPA 
1 k 



{U k ^,$) + $(d t -DV 2 )4>} 



(6) 



If the anomalous dimensions are not too large, the LPA 
already provides a good description of the effective po- 
tential as well as a rather accurate estimate of the expo- 
nent v governing the divergence of the correlation length. 
Since our main goal is to identify the non-perturbative 
fixed point governing the PC/GV transition in d = 1, we 
shall restrict ourselves, in what follows, to the LPA |18| . 

The NPRG equation for the effective potential, valid 
for all reaction-diffusion processes involving a single 
species, has been established in [19|. Studying a par- 
ticular model amounts to solving this equation in a sub- 
space defined by the symmetries of the problem, starting 
with the corresponding microscopic action 5*. The flow 
equation for the dimensionless potential u = k d+2 U k , ex- 
pressed in terms of the dimensionless fields t/j — > k~ d ijj 
and tp — > ip, reads (to lighten notations, we omit the 
implicit dependence of u on the running scale): 



d s u = -{d+2)u + dijju (lfi) - V d 



,(2,0),.(0,2) 



(l + u^ 1 )) 2 



where u <Jl ' p ^ — 



d n+p i 



and Vri 



-d+l. 



-d/2 



(7) 
In our 



d n ipd p if> ' a dT{d/2) 
problem, the effective potential must remain unchanged 
under the simultaneous transformations ip — > — ip and 
tp — » —ip, ( "parity-conservation" jZi symmetry of the 
PC/GV models) % H3- This leads to the existence 
of three quadratic invariant quantities, ip 2 , tp 2 , and iptp, 
from which all other invariant combinations of the fields 
can be built. Action JTJ can be expressed in terms of 
these invariants, but it also possesses additional features: 
the potential of the microscopic action S is proportional 
to 1 — ip 2 and vanishes for ip = (in the PC language). 
One can check that this structure is preserved by the 
renormalization flow, which further constrains the func- 
tional subspace in which the running potential evolves. 



FIG. 1: (Color online) Variation with d of the eigenvalues 
of the fixed points in the lowest-order LPA. Blue, dashed 
lines: pure annihilation fixed point Fa. Red, solid lines: non- 
perturbative fixed point F* whose eigenvalues are both nega- 
tive for | < d < 1.3784 . . . and complex-conjugated at larger 
d (only the real part is plotted). 



To summarize, the structure of the running potential 
defining our problem is 

u(ip, tp) = (l- ij} 2 )T{?p 2 ,ip$) with F(0, 0) = . (8) 

Postponing the numerical resolution of the partial dif- 
ferential equation J7J in the functional subspace defined 
above, we now perform a Taylor series expansion of the 
potential. In the absence of any information about the 
radius of convergence, the point around which the ex- 
pansion is performed matters, all the more so since we 
want to eventually truncate it. Here, we shall expand 
the potential around a non-negative solution of the sta- 
tionary equations of motion f^f = and ~ = 0. The 

former is always satisfied by ij> = 1 , while the solutions of 
the latter then correspond to either t/j — (the "origin" ) 
or to a (running) ip > ("the minimum"). We present 
results obtained around the origin, which involve lighter 
equations than those obtained around the minimum. 

The simplest truncation, that we now analyze in some 
detail, consists in keeping in u only the two coupling 
constants already present in S. Inserting this Ansatz in 
Eq.Q, we obtain the following non-trivial flows for the 
running const ants [2 lj: 



d s a 



-A (2 - d) + 2V d 



A 2 (l + 22g) 
(1 - 2a) 3 



-2<r + 6V d 



Act 



{\-2af 



(9) 
(10) 



This RG flow possesses three fixed points: the trivial, 
Gaussian, fixed point Fq (Aq = cr G — 0), the annihila- 
tion fixed point Fa (A a — f^r><?A = 0)> an d the non- 
trivial fixed point F* of coordinates: 

192 4 - 3d 



A* 



V d {28 - 3d) 2 



56 -6d 



(11) 



FIG. 2: (Color online) Flow diagram of the lowest-order LPA 
in d — 1 (as usual, arrows represent the RG evolution as s is 
decreased towards the "infra-red", macroscopic limit k — > 0). 



The Gaussian fixed point Fq, of eigenvalues (2, 2 — d) is 
relevant above d c = 2, where it coincides with Fa- For 
d € [§;2], Fa, whose eigenvalues are (d — 2, 3d — 4), is 
relevant. At d — |, Fa and the non-trivial fixed point 
F* coincide and exchange stability, so that F* is the 
relevant fixed point for d < | (Fig. Note that then 
a* > 0, and thus F* is in the physical region of parameter 
space, whereas it plays no role for the physics of reaction- 
diffusion systems when d > |. Note also that F* is not 
Gaussian in any dimension (at least at this order), and 
is thus out of the reach of any perturbative expansion. 
The PC/GV fixed point in d = 1 is thus F* , a genuinely 
non-perturbative fixed point. Its associated critical ex- 
ponent, given by the inverse of its negative eigenvalue is 
v = j= 12 - ~ 2.30. The flow diagram in this dimension 
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is shown in Fig. [3 The once unstable manifold of F*, 
connected to Fq, is the critical "surface" separating the 
absorbing and the active phases. The flow around Fa 
is rather peculiar: as d is decreased from |, the eigen- 
vector of Fa which is not parallel to the A axis rotates 
and becomes parallel, in d = 1, to this axis. Thus, Fa is 
degenerate for d = 1, since its two eigenvectors coincide 
(Fig. Q 22] . This implies in particular that every point 
flowing in the absorbing phase reaches Fa along the A 
axis. It is not clear to us, at this point, what might be 
the physical signature of this for microscopic models. 



TABLE I: Values of exponent v with the order n of the LPA 
truncation (see text). The column "min" refers to the minimal 
truncation (X,a). The last column is a conservative estimate 
taken from various Monte-Carlo simulations. 



n min 3 


4 


5 


6 


7 


8 9 


MC 


v 2.30 2.48 


2.20 


2.23 


2.11 


2.057 


2.015 2.0017 


1.85(10) 



We now report on the results obtained for truncations 
of the potential u that go far beyond the simplest trim- 
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cation described above. Of course, there exist many 
ways to organize a polynomial expansion of u around 
■0 = 0, ijj 2 = 1 in terms of the three quadratic invariants 
ipip, 1 — -i/> 2 , iji 2 which abides the Z 2 -symmetry of the 
PC/GV class. Equivalently, thanks to Eq.lJSJ, one can 
use any basis spanned by monomials in t/j 2 , ?/;■!/>, and we 
have tried several choices for it. In all cases, we stress 
that the qualitative picture unveiled at the minimal level 
is preserved at higher orders. It turns out that the fastest 
convergence of the exponent v is obtained with the ba- 
sis which is also the more transparent from the physi- 
cal point of view. At order "n" , we consider all pos- 
sible branching (n — 2m) A — > nA and annihilating re- 
actions nA — > (n — 2k)A involving at most n particles. 
These elementary reactions respectively correspond to all 
the terms (1 - 2m )(^)"~ 2m and (1 - $ 2k )^ n ^ n ~ 2k , 
and are anyhow ineluctably generated under renormal- 
ization. Table [I] shows the exponent v computed up to 
order 9. The convergence is rather good and we deduce 
that within the LPA v — 2 ± 0.1, already in fair agree- 
ment with the results obtained in numerical simulations. 
We finally report a rather surprising fact: at all orders of 
the truncation, the dimension at which F* and Fa cross 
is d = |. On the other hand we know, from the per- 
turbation expansion performed around d = 2 that this 
dimension is no longer | at two-loop order so that there 
is no reason to believe that this value is an exact result. 
Thus the deviation from | can only come from orders 
beyond the LPA in the derivative expansion. 

To summarize, we have shown, within the local poten- 
tial approximation of the NPRG, that the critical point 
of the PC/GV class of absorbing phase transitions is a 
genuinely non-perturbative fixed point, out of reach of 
perturbative methods. This result will have to be re- 
fined within higher-orders approximation of the deriva- 
tive expansion, so that a full set of critical indices can 
be estimated, not just the correlation length exponent v. 
Finally, ongoing work aims at a full numerical simulation 
of Eq. |0 (for the current problem and similar ones, such 
as the directed percolation class), which should allow to 
fully confirm these results and also offer access to other 
currently much debated reaction-diffusion problems, such 
as the pair contact process with diffusion [3| ■ 
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